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Abstract 

Requiring the two-Higgs-doublet model II to accommodate the 3a deviation in the muon anoma- 
lous magnetic moment imposes specific constraints on the Higgs spectrum. We analyze the com- 
bination of all the relevant, available, constraints on the model parameter space. The use of 
constraints from 6^57, the precision electroweak measurements of Rb, and the p parameter, 
together with exclusions from direct searches at LEP, give extremely severe restrictions on the 
model parameters. That is "almost enough" to kill the model altogether. The exclusion would be 
even stronger if the direct searches can be optimized to complement the other constraints, as will 
be discussed in details in this work. 

PACS numbers: 
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I. INTRODUCTION 



Physicists have been deUghted by the extraordinary success of the standard model (SM), 
while many got frustrated by the lack of experimental clues for the construction of the theory 
beyond. The most recent measurement on the muon anomalous magnetic moment [1] 
revealed a plausible deviation as large as Sex from the SM prediction, as suggested by the 
papers in Ref.[2]. While such a scenario is only a favorable, rather than unquestionably 
established, conclusion from the analyses, it has been taken by some physicists as a strong 
suggestion for physics beyond the SM. ^ Since many extensions of the SM are capable 
of giving rise to such a deviation, theorists, typically, would like to check the constraints 
imposed on the parameter space of a specific model. Such constraints are particularly 
interesting, because they are likely to give information not only on excluded regions but 
also on the predictions for where else should other evidences on the model be expected. For 
example, one expects to find a lower bound on the mass of a new particle playing the role on 
generating the extra contribution to a^. The information has especially strong implications 
on models that have a small number of parameters which are already stringently constrained 
by various precision electroweak data. We present here such a case study, illustrating how 
far such a 3o" deviation can take us. 

One of the simplest extensions of the SM is the two-Higgs-doublet model (2HDM) [3], 
which adds one Higgs doublet in addition to the one required in the SM. A generic 2HDM 
allows flavor-changing neutral currents (FCNC), which can be avoided by restricting the 
couplings of the doublets, say, by imposing an ad hoc discrete symmetry [4]. The most 
popular version, known as model II, has one Higgs doublet coupled to the down-type quarks 
(and charged leptons) and the second doublet to the up-type quarks. The physical content 
of the Higgs sector (assuming no CP violation) includes a pair of CP-even neutral Higgs 
bosons H and h, a CP-odd neutral boson A, and a pair of charged-Higgs bosons H^. The 
model fits in well with the criteria mentioned above for the a^^ result to have a very strong 
impact. We focus on the 2HDM II in this paper, answering the question of to what extent 
the model can survive a requirement of generating the 3 a deviation in a^, as suggested. 

The 2HDM II has been extensively studied in literature and tested experimentally. One 
of the most stringent tests is the radiative decay of B mesons, specifically, the inclusive decay 
rate of 6 ^ s 7, which has the least hadronic uncertainties. In the 2HDM, the rate of 6 ^ s 7 
can be enhanced substantially for large regions in the parameter space of the mass m^± of 
the charged-Higgs boson and tan/3(= tb/ti, where i\ and are the vacuum expectation 
values of the down- and up-sector Higgs doublets, respectively). An earlier analysis has 
already put a constraint on the charged-Higgs boson mass at m^± > 380 GeV[5] (see also 
Ref.[6, 7]). Updating the constraint while asking for the model to give rise to the deviation 



This 3c7 deviation was derived using the e+e data. If the r decay data are used, the deviation would be 
reduced to about 0.9cr [2], which, however, has model dependence and thus less reliable. 
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as suggested already imposes a strong and specific mass hierarchy between the pseudoscalar 
and the charged scalar. Electroweak precision data also have strong implications on the 
model. 

The 2HDM can explain the muon anomalous magnetic moment deviation with a light 
pseudoscalar boson A contributing via a 2-loop Barr-Zee-type diagram [8, 9]. Our interest 
here is in updating a previous analysis [9] and extending it to a comprehensive treatment of 
all the relevant constraints. The OPAL Collaboration [10] has recently published an update 
on their search for the Higgs bosons within the 2HDM framework. Since their result is more 
stringent than before, by combining the updated constraints from and other precision 
measurements, such as the p parameter, Ri,, and the b ^ sj rate, together with this new 
OPAL result and a study on Yukawa processes from DELPHI [11], we are able to limit the 
2HDM to a tiny window of parameter space with an, perhaps, uncomfortably large value of 
tan/3. All in all, we will piece together a story of the very stringently constrained 2HDM, 
almost to the extent of killing it altogether. We note that there have been previous analyses 
on the 2HDM using the electroweak precision data [12-14], to which we are partially in debt. 
Our study here can be considered an update, with a presentation along a different line. 

We are here talking about the intricate interplay of a few stringent constraints on the 
overall parameter space of the model. It is not a simple matter to illustrate results on the 
space of a large number (six, as discussed in the next section) of parameters. In our opin- 
ion, the best way to do it may depend on how each of the constraints really works. Our 
presentation of the constraints follows what we consider the most efficient way to appre- 
ciate the overall results. It may be not very conventional, but is considered particularly 
illustrative. Wherever explicit plots are shown, we are typically plotting two- parameter fits 
of one or more constraints, based on the usual x^-analysis. We will show < 4 regions, 
corresponding to 2a deviation limits, which we consider as "solution" regions — where the 
2HDM survives the particular constraints. We will also show regions where the model fits 
better than the SM, wherever appropriate. 

The organization is as follows. In the next section we briefiy describe the 2HDM (II) and 
the relevant parameters used in our analysis. In Sec. Ill, we look at the fo — > 57 rate and 
the — mixing, which require a heavy charged Higgs boson. In Sec. IV, we discuss the 
Higgs-sector contributions to and R^, and show the strongly complementary character of 
the two data and how that plays off in the 2HDM. We give our fits to the data, with the 
exclusions from DELPHI and OPAL further imposed. In Sec. V, we add the consideration 
of the p parameter. We conclude in Sec. VI. 

II. PARAMETER SPACE OF THE TWO HIGGS DOUBLET MODEL II 

We take the parameter space of the model as given by a set of six Higgs-sector parameters 

ruh , tUh , mA , mjj± , tan/3 , and a . 
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Here, the first four are masses of the physical Higgs states. The last one, a, is the real scalar 
Higgs mixing angle as defined in Ref.[3]. Wc would like to emphasize that we are taking 
this set of six parameters as mutually independent experimental parameters. From the 
theoretical point of view, one has parameters in the scalar potential from which the above 
can be derived. However, without supersymmetry or anything else to avoid the hierarchy 
problem, the masses for the Higgs states suffer from quadratic divergences. The tree-level 
relations among parameters are modified substantially by loop corrections which depend on 
the renormalization approach and the cut-off imposed. To stay away from such uncertainties, 
we do not discuss the scalar potential here, except noting that there are enough degrees of 
freedom in the model in general, and especially with the loop corrections taken into account, 
to allow us to take the above six parameters as mutually independent. Here we only consider 
the 2HDM-H without CP-violation, and are not interested in couplings among the Higgs 
states. It is then obvious that we do not have to consider more than the six parameters. 
It should be noted that any substantial CP violation in the Higgs sector that may largely 
invalidate our analysis here is ruled out by the electron electric dipole moment constraint. 

The four physical masses are direct experimentally measurable quantities. The other two 
parameters, tan/3 and the angle a, come into the game as effective couplings. The Yukawa 
couplings of h, H, and A to up- and down-type quarks are given by, with a common factor 
of —igmf/2Mw, 

tt bb T-T+ 

h: cosa/ sin/3 —sina/ cos/3 —sina/ cos/3 
H: sma/ sin/3 cosa/ cos/3 cosa/ cos/3 
A: —i cot/3 75 —i tan/3 75 —i tan/3 75 

while the charged Higgs H~ couples to t and b via 

btH" : [mt cot/3 (1 + 75) + tan/3 (1 - 75)] . 

2v2Mi4/ 

From the perspective of our study here, the couplings given above may be considered as 
defining implicitly the two parameters tan/3 and a. Other relevant couplings in our study 
are those to gauge bosons, as given by, 

, , , sin(/3 — a) 

hZZ : ig 714 ' 

coscV 

HZZ : ig M.^^?*^ 

COSt^VK 

, . ^ cos(/3 — a) , „ „ 

HAZ: -,i^^(,-,r 

2 cosOw 

H^H-Z : -ig^^ip-p'r, 
2 cos9w 

where p{h, H, and p'{A, H~) are the 4-momenta going into the vertex. 
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III. REQUIREMENT OF A HEAVY CHARGED HIGGS BOSON 



It has been well appreciated that B physics bars the 2HDM from admitting a relatively 
light charged-Higgs state. We first review the constraints here. The first important con- 
straint comes from the inclusive B — > Xg'-f result, and the second one comes from the B^ — B^ 
mixing. The essential point here is that without a direct source of FCNC, the charged Higgs 
mediates the only significant contributions to flavor-changing processes, in addition to the 
M^*-mediated SM process. The experimental data then allows us to bound the charged 
Higgs mass independent of the other Higgs states. 

The detail description of the effective Hamiltonian approach can be found in Refs. [15, 16] . 
Here we present the highhghts that are relevant to our discussions. The effective Hamiltonian 
for B — > Xg^ at a factorization scale of order 0{mb) is given by 



(1) 



The operators Qi can be found in Ref.[15], of which the Qi and Q2 are the current-current 
operators and — Qq are QCD penguin operators. (^77 and Q^g ^-re, respectively, the 
magnetic penguin operators specific for b ^ s'j and b ^ s g. Here we also neglect the mass 
of the external strange quark compared to the external bottom-quark mass. There have 
been more recent analyses [7, 17, 18] on 6 ^ 57 involving the NLO and other corrections, 
but the LO treatment here is sufficient for our purpose to put a lower bound on the charged 
Higgs mass mH±, which is then used in the central part of our analysis. 

The decay rate of i? — > Xg'y normalized to the experimental semileptonic decay rate is 
given by 

where f{z) — 1 — Sz'^ + 8z^ — z^ — 2Az^\nz. The Wilson coefficient C-j^{mij) is given by 

C,^{^J) = ?7^C'77(M^) + I (?7^ - ?7^) C^g{M^) + C2{M^)j2hiV''' , (3) 

1=1 

where rj — as{Mw)/ois{ij). The Oj's and /ij's can be found in Ref. [15]. The coefficients 
Ci{Mw) at the leading order in 2HDM II are given by 

C,(M^) = (j = 1,3,4,5,6), (4) 

C2(M^) = 1 , (5) 
C,,{M^) = - ^ coi'p - B{y,) , (6) 

Csg{M^) = - ^ cot^/3 - E{y,) , (7) 
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where Xt — rn^/M^, and yt — m1/m'^±. The Inami-Lim functions [19] are given by 



A{x) 
B{y) 

D{x) 

E{y) 



X 



12(a;-l)3 2{x-lY 
5y-3 (3y-2)lny 



I2{y - 1)2 6(y - 1)3 

^x"^ — 5x — 2 3a; In a; 



X 



= y 



4(a;-l)3 2(a;-l)^ 



+ 



\ny 



(8) 
(9) 
(10) 

(11) 



U{y-iy 2{y-lf. 
The most recent experimental data on 6 — > 57 rate has been reported [20], giving 

B{b S7)|exp = 3.88 ± 0.36(stat) ± 0.37(sys);°:^^(theory) . 

The most updated SM prediction is [21] 

5(6 ^ s 7) Ism - (3.64 ± 0.31) x 10"^ , 

which agrees very well the data. Both the experimental data and the SM prediction have 
been extrapolated to the total branching ratio. Therefore, there is only a httle room for new 
physics contributions. The constraint on new physics contribution is, explicitly. 



AB{b -^s^) = B{h^ S7)|exp - Bib ^ S7)|SM = (0.24;°j^) x 10 



(12) 



where we have added the various errors of the experimental data in quadrature. (Note that 
the theory error quoted in the experimental data is larger than the one quoted by the SM 
prediction. We take the more conservative value.) 

The quantity that parameterizes the — B^ mixing is 



where[22] 



Xd = 



F ITA* |2|TA |2 f2 



67r2 



(13) 



1 + 



3 — 9a; Gx^ log a; 



ix-1) 



xy cot^P 



(x-l) 



{4:Z - 1) logy 



3 logx 



+ 



X - 4 



2(1 - y)2(l - z) 2(1 - a;)2(l - z) 2(1 - a:)(l - y) 



xy cof^P \ 1 + y 2y log y 
(1^+ (1^ 



with X = m1/M^,, y = m^/m|±, z = M^^,/m^±, and the running top mass rrit — mt{mt) — 
166 ± 5 GeV. The experimental value is [23] 



Xd = 0.755 ± 0.015 . 



(14) 
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We use the following input parameters [23]: \VtbV*a\ = 0.0079 ± 0.0015, f^Bj, = (198 ± 
30 GeV)2(1.30 ± 0.12), = 5279.3 ± 0.7 MeV, r/^ = 0.55, and = 1.542 ± 0.016 ps. Note 
that the value of jV^bV^^I is in fact determined by the measurement of x^- Now we can use 
the data to constrain the new contribution from the charged-Higgs boson. 

The two constraints discussed above are quite stringent, giving a lower bound on m„± 
close to 500 GeV at 95% C.L. for intermediate and large values of tan/9. The result is 
illustrated in Fig. 1. The heavy charged Higgs mass means that the state pretty much 
decouples, playing a little role in the contributions of the 2HDM to quantities like and 
Rb, which we turn to in the next section. 

IV. VS Rb 

The most recent data on the indicates [2] (see also the footnote # 1) 

Aa^ = a^^P - af^ = (33.9 ± 11.2) x 10"^° . (15) 

The result shows a 3cr deviation to be explained by new physics. Adopting the view that 
the problem is real and demands new physics contributions, we will see that it has 
a strong and definite implication on the Higgs spectrum of the 2HDM. In fact, we had 
performed an analysis [9] along the line for the earlier data. The major point is that a 
light pseudoscalar, together with a large tan/3 value, is required to explain the positive Aa^ 
contribution via a two-loop Barr-Zee diagram [24]. A real scalar contributes in the negative 
direction. To avoid a cancellation, the real scalar mass has to be heavy. We will show here 
that such a mass splitting is strongly disfavored by the allowed contribution to Rb, for which 
the experimental data agrees well with the SM prediction. The two constraints are hence 
strongly complementary. 

It has been emphasized in Ref. [8, 9] that for the Higgs boson mass larger than about 
3 GeV, the dominant Higgs contributions to actually come from the two-loop Barr-Zee 
diagram with a heavy fermion (/) running in the upper loop. A m^/m^ factor could easily 
overcome the a /An loop factor. In our calculation here, we include all one-loop contributions 
and all two-loop Barr-Zee-type contributions with an internal photon and a third-family 
fermion running in the loop. The latter diagrams with the bottom and tau loops are strongly 
enhanced by tan/9. If the internal photon was replaced by a or a. Z^, the contributions 
will be much suppressed (see [25], for examples). The case is in particular strongly 
suppressed, partly as a result of the fact that the Higgs boson has to be H^, the mass of 
which we have shown above to be heavier than 500 GeV. The only other important diagrams 
of the Barr-Zee type are the SM diagrams, such as the one with the replacing the heavy 
fermion. We neglect the small "extra" contributions from such diagrams[26] because of the 
small difference between the Higgs boson mass used here and that used in the Ref. [27]. 
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Explicitly, we first write the fermion couplings of a neutral Higgs mass eigenstate 0° as 

C^^f = -A, + Af , (16) 



where A/ ^ and Af ^ are the effective scalar and pseudoscalar couplings explicitly given 
in Sec. II, and v — 246 GeV. The two-loop photon Barr-Zee diagram contribution from 0°, 
with a heavy fermion / running in the second loop, is given by 



(17) 



where 



2 Jo x{l — x) — z z 



/ X 1 /"^ , 1 . x(l- x) 

g{z) = -z dx — ^ In ^ ^ ; 18 

2 Jo x{l — X) — z z 

NJ represents the number of color degrees of freedom in /, and Qf its electric charge. Here, 
we have three scalars and no CP violating mixing is assumed. The real scalars h and H give 
negative contributions only (from the second part) while a pseudoscalar A gives a positive 
contribution only (from the first part). The diagrams with the b and r loops are tan/3 
enhanced. 

We plot in Fig. 2 the 2a range of solution to on the plane of the pseudoscalar mass 
ijIa versus tan/5, considering only the pseudoscalar contribution. Note that while the re- 
gion below the solution band is excluded, the solution in the region above the band may 
be admissible when the is compensated by some negative contributions from the real 
scalar (s). We also superimpose on the plot the excluded region from the DELPHI study on 
Higgs Yukawa processes in the 46 and 262r final states [11]. We can see that one obtains 
lower bounds on uia and tan/3 as 26 GeV and 30 respectively. ^ We had given in Ref. [9] 
plots of the solution regions on the ruh-mA plane for the old data with specific values of 
tan/3 and the scalar mixing angle a. We will present similar results here with, however, the 
complementary constraint included. We will see that not much area of the parameter 
space can survive the combination of both and data. 

The current it!;, measurement is given by [29] 

i^fP = 0.21646 ± 0.00065 . 



^ The plot is similar to the one given in Ref. [14], in which some more constraints are superimposed. We 
include here only the important ones. Note that the Tevatron exclusion region claimed in the paper 
is not used here. The exclusion result was from Ref. [28], which is an analysis based on the minimal 
supcrsymmetric standard model. The result should not bo directly applicable to the present case of 
2HDM. We do not find a similar study on the Tevatron data based on the 2HDM. 
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With Rf^ = 0.215768, we have 

ARb = - Rf^ = 0.000692 ± 0.00065 . (19) 

The ARh contributions in the 2HDM are given, for example, by formulas in Ref.[30]. The 
charged Higgs contribution is always negative. On the other hand, there is a window of 
parameter space for the neutral Higgs contributions to be positive. From our discussions 
above, we need a light pseudoscalar and have to live with a heavy charged Higgs. We are 
therefore more interested in the neutral Higgs contributions. 

Let us first focus on the contributions from the lighter Higgs bosons, the pseudoscalar 
A and the real scalar h, pushing the other scalar H to the heavy-mass limit together with 
H^. The scenario will actually be well justified by our discussion on the constraint of the p 
parameter in the next section. The smallness of AR^ contributions admitted here generally 
disfavors a large mass splitting between nih and rriA- This is in contrast to the requirement of 
a positive contribution to a^. Since the SM result (with Mhsm 115 GeV) now represents 
a 3(7 deviation in a^, better fits to the combined a^-R}, data are possible from the 2HDM. 
We illustrate some such fits in Fig. 3. In the figure, we take the case of tan/? = 58 and 
check various values of the Higgs mixing angle a. Here, and in the discussion below unless 
specifically stated otherwise, we stick to r?7jf± = 500 GeV and = 1 TeV. The exact value 
of itlh does not matter at all here, one, however, should note that the charged Higgs boson 
still gives a contribution of —2.32 x 10^^ to which is about in strength. Bearing this 
in mind, it is easy to estimate from our plots the slight shift in each of the admissible region 
as m^± is being pushed towards the decoupling limit. 

Each of the plots in Fig. 3 gives ±2cr limits for and R}, fits, with darker shaded regions 
indicating the solution of interest defined by a total of 4 or less. Also marked in the plots 
are regions with a total less than the SM value of 10.3 (the sum of and i?^). The 
purpose of showing the area with a total < x^(SM) is to indicate the region of parameter 
space that can fit the and Ri, better than the SM, other than the decouphng hmits. From 
now on, we concentrate on the dark area of a total x^ < 4 as a valid solution to the and 
Rb data. We can see that there are no solutions for — | < a < |. While a larger magnitude 
of the a looks more favorable, the best solution, with higher Higgs masses (m^ especially), 
stay close to | sina;| = 1, inclining more towards the negative sign. The most favorable range 
is around — | < a < — ^ [c/. plots (c) and (d)]. The plots in Fig. 3 illustrate well the trend 
of the changes in the a^-Rh solution regions with variations in a, which is quite generic for 
tan/3 around and larger than 50. A higher rrih solution is preferred as the light mass solutions 
are easily killed by searches at LEP, including the DELPHI exclusion used in Fig. 2 and a 
particularly focused 2HDM analysis from OPAL [10], to be discussed below. In fact, as we 
will illustrate below, if any part of the solutions to the a^-R^ fits survive the exclusions from 
LEP, it is more or less the part with a high enough nih value. ^ 



^ In these plots, we used the standard Rh formulas for 2HDM as available in the reference [30] quoted. The 
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For instance, the DELPHI exclusion we used in Fig. 2 obviously kills quite a part of the 
above solutions. At tan/5 = 58, the lower bound on is actually about 40 GeV. It is 
particularly interesting to check the case of a relatively small tan/5. We show in Fig. 4 the 
case of tan/5 = 40, in which the DELPHI almost kills all solutions. In the figure, only a tiny 
window survives at sina = This is, however, at a nih (as well as m^) value too small to 
survive the OPAL exclusion discussed below. From the same study by DELPHI, the Yukawa 
processes were also used to impose bounds on the m/i[ll]. Note that the plots in the Ref. [11] 
give bounds on Higgs masses versus the b — r coupling enhancement factor, which is simply 
ta,n(3^ B(A bb, t+t~) for the pseudoscalar case, but has an added | sinal dependence for 
the case of the scalar h. In the range of uia values of interest, however, the ruh bound is 
typically superseded by the OPAL exclusion. Note that the solution with nih at a few 
GeV is also inadmissible, from the consideration of T-decay[9] and otherwise. Figure 4 also 
illustrates that a positive value of a, around tends to give a^-Rb solutions with the largest 
niA- This is mainly a result of the constraint. At a fixed m/j, the contribution of the 
scalar through the 2- loop Barr-Zee diagram to is suppressed by | sinal, and thus allows a 
larger m^. A slight asymmetry in the cases of positive and negative a values comes in as a 
result of the different way the 1-loop contributions go. 

We should also point out that we find no solutions to the a^^Rb fit at all for much lower 
tan/?. The solutions shown in Fig. 2 simply produce a mh-rriA splitting too large to 
accommodate the i?5 constraint. Therefore, the solutions to the a^-Ri fits start to emerge 
as tan/5 approaches a large enough value, not much below 40. To get solutions with high 
enough masses for the Higgs bosons so as to survive the exclusion limits from the LEP 
searches, one will have to get to a higher and higher tan/5 value. To check the details, we 
first turn to the powerful exclusions from OPAL. 

Based also on the LEP data, OPAL has been pubhshing Higgs-search analyses specifically 
focused on the 2HDM. Here, we use the most recent results available [10, 31]. The results 
exclude a region at the lower left corner of the m^i-m^ plane for each specific value of the 
mixing angle a (four explicitly shown). As presented in Ref. [10], however, the results are 
not tan/5 specific. We show in Fig. 5 the solution regions from a^-Rb with the OPAL and 
the above mentioned DELPHI excluded regions superimposed. Here, we show the cases 
of tan/5 = 50 and 58 for two a values, — | and ±|. The latter are chosen as they are 
among the a values for which the OPAL paper gives the explicit exclusion. The exclusion 
for non-specific a values only presents a substantially weakened result to be of interest here. 



formulas have not taken into consideration the possible tree-level decays of the boson into a pair of 
Higgs states. Such very light Higgs regions are typically not of interest to theorists, as they would be 
easily excluded by direct searches. By the same token, we use the formulas as they are and show the 
a^-Rb fits without bothering about the tree-level modifications of Ri, from such Higgs channels, leaving 
such fake solutions to be taken care of by the exclusions from LEP searches to be discussed below. We 
thank A. Sanda for alerting us to explicitly address the issue here. 
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The two a values are also close to or within the range of value for an optimal a^-Rb fit, as 
illustrated in Fig. 3 above. 

In fact, the OPAL exclusions are based on the tan/5 value in the range 1 — 58, given only 
at four values of a, which we adopted here for the plots. An cxchidcd point is one excluded 
at all value of tan/? within the range. The excluded ranges may hence be extended at each 
specific value of tan/5, with more detailed analysis of the data[32]. In fact, one would expect 
the enhanced couplings at a larger tan/3 generally push the exclusion regions towards higher 
masses. One can see in the plots that the a^-Rb solution regions are largely cut off by the 
presented OPAL exclusions in general. Actually, nothing survives in the illustrated plots 
for tan/3 = 50 [c/. plots (a) and (b) of Fig. 5]; and it is also quite obvious that the same is 
true for the case of tan/3 = 40. Figures 4 and 5 together clearly illustrate the general trend 
of how the increase in tan/3 gives better results. While it does not do much for the case of 
q; = — ^, atQ; = ±|, the admissible ruh and ruA values are pushed to be high enough to 
escape the OPAL exclusion at tan/3 > 50 [c/. plot (c) of Fig. 5]. However, we are not bold 
enough to say for sure if no solution survives at tan/3 = 50 though. The optimal case giving 
the largest rrih value for the a^-Rb is likely to be around a — — at which the present 
result of OPAL did not show its greatest strength. 

We pick tan/3 = 58 for the above detailed illustration of the a^-Rb fits because it gives the 
most favorable case within the limit of a direct application of the OPAL exclusion. There is 
apparently a surviving window in the parameter space. Combining our present study with 
a refined version of the OPAL analysis to focus on our a^^-Rh solution regions, especially 
the upper blots as shown in plots (a) and (c) of Fig. 5, will certainly be very interesting. 
Exclusions have to be checked for each specific value of tan/3 and that of a. Such a study 
will further narrow down the survival parameter space regions of the 2HDM, especially with 
further improved exclusion limits. We are told that the LEP data actually allows such an 
improvement [32]. A very tantalizing question is if the current constraints are actually strong 
enough to kill the model altogether! 

The OPAL analysis is hmited to tan/3 value at or below 58, while the DELPHI analysis 
stops at 100. The very large tan/3 region is theoretically unfavorable and may provide 
practical problems due to the much enhanced 6-quark Yukawa coupling, which signals a 
breakdown of the perturbative treatment. Nevertheless, we will include some results from 
such uncomfortably large tan/3 values, and urge the OPAL group to push on a bit further 
in their analysis. 

For tan/3 > 58, we again illustrate some results for a = — ^ and ±| in Fig. 6, in which 
we still put in the (no longer exactly valid) OPAL exclusion from the tan/3 = 1 — 58 scan for 
reference. For the a = ^ case, the a^-Rb solution regions never rise above rUh = 60 GeV 
(as also for the smaller tan/3 cases), and partially excluded by DELPHI. There seem to be 
good enough reasons to believe that such low Higgs mass regions should be excluded by the 
available LEP data if an analysis along the OPAL hue is performed. The a = ±| case is 
better. The surviving region seen at tan/3 = 58 moved further to the right, towards larger 
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rriA, and further up a bit for larger rrih- The rise in nih actually more or less saturates at 
tan/3 = 100, and falls for even larger values of tan/?. Such a region still lives at the boundary 
of the DELPHI exclusion, but would still have a surviving part even if the OPAL exclusion 
can still be imposed. The region is, however, shrinking with increases in tan/3, due to a more 
fine-tuned a^-i?;, solution. One should also bear in mind that the (upward) shifting, and 
hence slight enlargement, of this solution region with a further increase in the m^± value. 
For this purpose, we give the charged-Higgs contribution to Rj, in Fig. 7. 

Let us summarize our results so far. We have seen that combining the suggested require- 
ment of producing a definite positive contribution to while keeping a limited deviation 
from the SM Rb result is an extremely stringent constraint on the parameter space of the 
2HDM. When the available direct experimental search results from the LEP experiments 
are further implemented, there is at most a tiny window of parameter space that can sur- 
vive. The apparently surviving region from the above discussions is restricted to very large 
values of tan/3. In fact, it may be already uncomfortably large, inviting the problem of the 
perturbativity of the Yukawa coupling of the b quark. As for the mixing angle a, it is being 
pushed close to the — | < a < — ^ region. All these are based on a strong mass splitting 
between the charged Higgs and the light scalars, with the pseudoscalar A lying typically 
below 80 GeV and the scalar h below 140 GeV. 

Recall that we have essentially decoupled the heavy Higgs boson H in the above analysis. 
We promised to justify this as a physics requirement. We first noted that a not too heavy H 
would more or less add to the effect of the other real scalar h in the contributions to and 
Rb. So, we expect it to ask for a larger mass splitting when fitting is concerned, but a 
smaller mass splitting to fit Rb- In another word, further tightening of the apparent solution 
window. In the section below, we will show that fitting another precision EW parameter, 
the p parameter, actually does require sending mn to a very large value, indeed well beyond 



V. THE p PARAMETER CONSTRAINT 

The parameter p was introduced to measure the relation between the masses of and 
Z° bosons. In the SM p = M^/M^ cos^^ = 1 at tree- level. However, the p parameter 
receives contributions from the SM corrections and from new physics. The deviation from 
the SM prediction is usually described by the parameter po defined by [33] 



where the p in the denominator absorbs all the SM corrections, including the corrections 
from the top quark and the SM Higgs boson. By definition, po = 1 in the SM. Sources of 
new physics that contribute to po can be written as 



Po = 




(20) 



pM|cos26U ' 



Po = 1 + Ap; 







new 



(21) 
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where Ap"^"^ = Ap^"^'^ - Ap^^-^'^^ in our case. Note that since the two-doublet Higgs sector 
(in the 2HDM) is employed here to replace the SM Higgs, the latter contribution to Ap has 
to be subtracted out. 

The most recent reported value of po is [23] 

Po = 1.0004 ± 0.0006, (with 714^^ fixed at 115 GeV) . (22) 

In terms of new physics the constraint becomes: 

Ap^*^™ = 0.0004 ± 0.0006 . (23) 

In 2HDM Ap receives contributions from all Higgs bosons given by [3, 12] 

Qj r 

^'^^ = , • 2T ^("^^' "^+) + - [F{m„+,mh) - F{m^, m^)] 

47rsm OwJ^w 

+ sin^/? - a) [F(mff+, mn) - F{mA, nifi)] 

+ cos\p - a)Ap'^\m^) + sm\p - a)Ap'^\mh) , (24) 



where 



F{x,y) = -X +-y - l^^TT^^^'Si-^ j =F{y,x), 



47r sm^0„M3, 



3F{M, M^) - 3F{M, M,) + i(M| - M^) 



(25) 



Let us take a closer look into the implication of the formulas above. First of all, we note 
that Ap^'^M) has a negative value with magnitude increasing with M. As to be expected 
from above, the value is about —0.0004 at M = 115 GeV. It has a relatively mild variation, 
and does not go beyond —0.005 even as M gets to lOTeV. In direct contrast, the other 
contributions to Ap^^""^ in the above formula are very sensitive to the masses involved. The 
F{x,y) function is always positive, vanishes only at x = y, and increases with a faster and 
faster rate with the splitting between x and y. In a typical scenario that is of interest here, 
we expect the pseudoscalar to be the lightest Higgs state with a quite heavy charged Higgs. 
That makes the contribution from the first term [involving F{mA, Tn„+)] large; indeed of order 
0.01 for satisfying the lower bound from b ^ S'~f and B — B mixing. To get the required 
almost zero value of Ap^^°^, we need some negative contributions from the terms involving 
[F(m^+,m/i) — F(mA,m/i)] and [F(m^+,'m«) — F(m4,m^)]. And the solution is obviously 
a fine-tuned one. Consider a case of degenerated Higgs real scalars, tUh = ruh = M. The 
{(3 — a) dependence in the formula is removed as the sine-square and cosine-square parts are 
combined. We need [F{mH+,M) — F{mA, M)] to have a negative value close in magnitude to 
that of F{'mA, Tau+). The former is obviously positive roughly when M is closer to 7714 than 
771^+, and negative when it is the other way round; and it can be larger than F(mA,r?7jy+) 
only for M > m^+. The {(5 — a) dependence comes back in the generic situation, with a 
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mass splitting between the two real scalars. Obviously, at least one of them, by definition 
H, has to be heavy. However, they cannot be both heavy, because the Rb constraint does 
not allow only a light pscudoscalar giving a substantial contribution. Hence, this additional 
requirement of a limited splitting between hIa and nih clearly suggests a large m^, typically 
larger than m^+. The larger the mjj value, the smaller the sm\P — a) is required. Admissible 
solutions, though fine-tuned, can be obtained so long as the required sm\(3 — a) falls into 
the legitimate interval. A large splitting between and nih also resulted in a very narrow 
range of admissible (/3 — a) values, and hence the a values at a fixed tan/3 of interest. 

Our numerical results corroborate well with the above analytical discussions. We illus- 
trative our discussion with a plot in Fig. 8. Here, we take a "surviving" solution point to 
the a/^-Rb fits and perform a further fitting together with the p parameter by varying mn 
and m^+. The extremely fine-tuned nature of the solution is well-illustrated by the very 
narrow bands. Moreover, ttIh is always more than twice of for an | sinal larger than 
0.8. Note that sinct of —0.8 and —0.92 roughly correspond to an a value of and 
respectively. The basic features remain if other solution points are taken, hence, we refrain 
from showing more plots. 

VI. CONCLUSIONS 

The requirement for the 2HDM 11 to give rise to the suggested +3a deviation in the muon 
anomalous magnetic moment demands a light pseudoscalar, preferably around 40 GeV for 
tan/3 < 58. The FCNC constraints, in particular the 6 — > S7, require a heavy charged Higgs 
beyond 400 GcV. This large mass splitting in the Higgs mass spectrum is difficult to be 
accommodated by the precision EW measurements. In particular, the p parameter constraint 
then admits only very fine-tuned solutions, with cancellation from opposite contributions 
good to one in a hundred, favoring heavy real scalars. The fine-tuned nature of the solutions, 
while making many physicists uncomfortable with the model, is not in itself a good enough 
reason to pronounce the death of the model. If one has reasons to be confident about the 
correctness of the model, one would say that the available constraints are just strong enough 
to pin down for us the values of the unknown model parameter. A good example of such a 
situation is given by the pinning down of the top mass value from precision EW data prior 
to the experimental discovery of the top quark. 

Nevertheless, we have not had much of a reason to believe in the correctness of the 
2HDM. In fact, in our analysis here, we focus more on the simultaneous fits to Aa^ and 
Rb- Having only the pseudoscalar contribution dominated the corrections to Rb is fatal. 
Hence, we require a relatively light rrih, while pushing to way beyond in order to 
satisfy the p parameter. Even then, not much of a solution to the a^-Rb fits survives the 
direct search exclusion. Here, we require a total of 4 or smaller for the a^-Rb fits to claim 
a good solution. For roughly tan/3 < 40, no solution to a^-Rb fits survives the DELPHI 
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exclusion. In fact, there is no solution to the a^-Rb fit for a tan/3 value quite a bit smaller 
than 40. For tan/3 value from around 40 to a bit beyond 50, solutions surviving the DELPHI 
exclusions exist, but only to be killed by the OPAL exclusions. For even larger tan/5, the 
a^-i?5 solutions surviving both the DELPHI and OPAL exclusions started to emerge. This 
is very much restricted to an a value in the range — | < a < — Solution regions shrink 
fast outside the range as the ruh values given by the a/^-Rb fits drop towards OPAL exclusion 
bound. 

In summary, under the strong restriction of the available constraints, we show that only 
a very tiny window of apparent solutions exist close to the limit of tan/3 < 58. If the OPAL 
group could tailor their analysis of the LEP data to focus on the apparent solution window 
as shown here, we would be able to have a more definite conclusion. It looks to us that 
the solution window will be shut down quite substantially. So, we have "almost enough" 
constraints to kill the model altogether. For tan/3 beyond the 58 limit, our hands are tied 
at the moment by the unavailability of the strong LEP exclusion results as presented by 
OPAL. The very large tan/3 values certainly make many of us uncomfortable though. It is 
theoretical undesirable as limited by the blowing up of the bottom Yukawa couplings. 
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FIG. 2: The 2cr allowed region in the (mA, tan/3) plane due to the data. Here only the 
pseudoscalar contribution is considered. The DELPHI excluded region is represented by the thick 
line. 
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FIG. 3: The 2a allowed regions in the (m^, ruh) plane due to the constraints of and Rb for 
tan/3 = 58. The smaller dark region is where the total is less than 4 while the lighter region is 
where the total is less than the x^(SM) = 10.3. 
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(a) tan|3=40, a=37c/8 



(b) tanp=40, a=±jt/2 




FIG. 4: The 2a allowed regions in the (mA, nifi) plane due to the constraints of and Rf, for 
tan/3 = 40. The smaller dark region is where the total is less than 4 while the lighter region is 
where the total is less than the x^(SM) = 10.3. Here the lower mass limits on and ruh from 
DELPHI are shown. 
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FIG. 5: The 2o" allowed regions in the (m^, rrih) plane due to the constraints of and for 
tan/3 = 50 (a,b) and for tan/3 = 58 (c,d). The smaller dark region is where the total is less than 
4 while the lighter region is where the total is less than the x^(SM) = 10.3. Here the lower mass 
limits on and m/j from DELPHI, and the excluded region from OPAL are shown. 
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FIG. 6: Same as Fig. 5, but for tan/3 = 80 (a,b) and for tan^ = 100 (c,d). Note that the OPAL 
exclusion regions put in is actuahy for tan/5 < 58, hence no longer directly applicable here. In lack 
of the corresponding apphcable results, they are kept here for reference. 
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FIG. 7: The charged-Higgs contribution to in units of 10 ^ for tan/3 = 40 — 100 (from the 
to the bottom). 
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FIG. 8: The allowed region in the plane of (m^,m^+) due to constraints of a^, Rb, and the p 
parameter. Here the = 40 GeV and mh = 100 GeV for tan/3 = 58 are chosen from the allowed 
darker region in Fig. 5(c). 
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